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ABSTRACT 

Intracluster magnetic fields with ~ yuG strength induce Faraday rotation on 
the cosmic microwave background (CMB) polarization. Measurements of this 
effect can potentially probe the detailed structure of intracluster magnetic fields 
across clusters, since the CMB polarization is a continuously varying field on the 
sky, in contrast to the conventional method restricted by the limited number of 
radio sources behind or inside a cluster. We here construct a method for extract- 
ing information on magnetic fields from measurements of the effect, combined 
with possible observations of the Sunyaev-Zel'dovich effect and X-ray emission 
for the same cluster which are needed to reconstruct the electron density fields. 
Employing the high-resolution magneto-hydrodynamic simulations performed by 
Dolag, Bartelmann & Lesch (1999) as a realistic model of magnetized intraclus- 
ter gas distribution, we demonstrate how our reconstruction technique can well 
reproduce the magnetic fields, i.e., the spherically averaged radial profiles of the 
field strength and the coherence length. 



1. Introduction 

The origin and evolution of cosmic magnetic fields are still unclear and outstanding 
problems. Various observational techniques have consistently revealed that most clusters 
of galaxies are pervaded by magnetic fields of B ~ 0(l)/iG strength (see e.g., Carilli & 
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Taylor 2002 for a review). Recently, Clarke, Kronberg & Bohringer (2001) have again drawn 
a robust conclusion that the intracluster hot plasma typically has a B ~ 5 — 10 /zG field 
assuming a coherence length of lOkpc from the Faraday rotation measurements of 16 low-z 
[z < 0.1) clusters selected to be free of unusual radio halos. Since the rotation measure arises 
from the integration of the product of the electron density and the line-of-sight component 
of magnetic fields, it is generally difficult to extract information on the magnetic fields only 
from the rotation measure without introducing any assumptions on the gas distribution and 
the field configuration. 

The first systematic study of the structure of magnetic fields over a single cluster was 
performed by Kim et al. (1990) based on the rotation measures of 18 sources close in 
angular position to the Coma cluster, giving the estimation of B ~ 2yuG with the coherence 
length of lOkpc which is indicated from the magnetic field reversal scale. In a subsequent 
study, Feretti et al. (1995) discovered smaller coherence lengths down to 1 kpc from the 
rotation measures of the extended radio galaxy near the Coma cluster center, whereby the 
field strength estimation was modified to B ~ 8/zG to explain the measured rotation angle. 
Thus, it is crucial for estimating the magnetic field strength robustly to determine the 
magnetic field coherence length which does not necessarily match the coherence length of 
the rotation measures. To study the detailed structure of the intracluster magnetic fields 
for any clusters, the high-resolution measurements of the Faraday rotation should at least 
be performed. However, there are some limitations for the conventional methods because of 
the lack of the number of radio sources behind or inside a cluster and possible contributions 
of the intrinsic Faraday rotation. Moreover, we should break the degeneracy of the rotation 
measure between the electron density and the magnetic field strength using some additional 
information for estimating the magnetic field strength. 

Recently, Takada, Ohno & Sugiyama (2001) proposed that the magnetized intracluster 
gas similarly induces a Faraday rotation effect on the linearly polarized radiation of the 
cosmic microwave background (CMB) generated at the decoupling epoch of z ~ 1000 (see 
e.g., Hu & White 1997 for a study of the primary CMB polarization). Hereafter, to dis- 
tinguish the origin of the CMB anisotropies, we use "primary" to indicate the anisotropics 
from the recombination epoch while "secondary" to those from the cluster formation epoch. 
They also calculated the angular power spectra of this secondarily induced polarization un- 
der the simple assumption of a uniform field strength with ~ O.lyU G across a cluster and 
suggested that the measurements could be used to set constraints on the average properties 
of the intracluster magnetic fields. Cooray, Melchiorri & Silk (2002) computed the secondary 
power spectrum including a circularly polarized contribution characterized by the Stokes-V 
parameter that is induced by possible relativistic plasma in clusters via the Faraday rotation. 
However, as a more interesting and realistic possibility, one may imagine that the secondary 
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effect can be in principle used to map the detailed structure of the magnetic fields in an 
individual cluster, since the CMB polarization is a continuously varying field on the sky. In 
this paper, therefore, we study a method for reconstructing the magnetic fields from mea- 
surements of the Faraday rotation effect on the CMB polarization, combined with accessible 
observations of the Sunyaev-Zel'dovich (SZ) effect and X-ray emission for the same cluster 
to reconstruct the electron density distribution. For this purpose, it is crucial to consider 
a realistic magnetic field configuration as well as a plausible gas distribution as suggested 
by the currently favored formation scenario of galaxy clusters. Hence, as for a model of 
the magnetized cluster, we employ here high-resolution magneto-hydrodynamic simulation 
results performed by Dolag, Bartelmann & Lesch (1999), whereby we can simulate maps of 
the CMB polarization including the Faraday rotation effect as well as maps of the SZ effect 
and the thermal X-ray emission. Using those simulated maps, we demonstrate how well our 
method can reconstruct the magnetic fields. In particular, we try to clarify the relationship 
between the coherence lengths of the magnetic fields and the rotation measure. The coher- 
ence length should give a new insight into the nature and evolutionary history of magnetic 
fields. For example, if the magnetic fields have the coherence scale as large as or larger than 
the cluster size, the seed fields should be generated at the early stage of the universe (see 
e.g., Grasso and Rubinstein 2001 for a review), while smaller coherence lengths may imply 
the seed fields originated from galaxies within the cluster (e.g., Kronberg 1996). 

This paper is organized as follows. In §2, we refer to the cluster models which are used 
to demonstrate the magnetic field reconstruction. In §3, we briefly review the primary CMB 
polarization map generated at the decoupling epoch and the Faraday rotation effect. In §4, 
assuming spherical symmetry of clusters, we show a method for reconstructing the density 
fields of the clusters from SZ effect and thermal X-ray emission which are directly calculated 
from simulated clusters. In §5, we develop a way to reconstruct the coherence length and 
strength of the magnetic fields combined with the previously reconstructed electron density 
fields. Finally, §6 is devoted to a summary and discussion. 



2. Cluster models 

We employ high-resolution magneto-hydrodynamic simulations performed by Dolag et 
al. (1999) as a realistic model of magnetized clusters (see also Dolag et al. 2001a,b; 2002). 
The simulations start from the seed magnetic field with 10 _9 G strength at redshift z = 20, 
which corresponds an upper limit to IGM fields set by Faraday rotation measurements of 
high-z radio loud QSOs (Kronberg 1996; Blasi, Buries, & Olinto 1999). They considered the 
homogeneous or chaotic field configurations for the initial seed fields, which are motivated 
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Fig. 1. — Average electron density along the line of sight, (n e (l/cm 3 )) = l/(2L) J L L dsn e 
for each cluster (A, B, and C). The axes x and y show the angular scale in degrees. 



Table 1. Cluster Models 



Cluster A ClusterB ClusterC 



z 0.125 0.125 0.452 

r vir [Mpc] a 2.7 2.7 2.1 

M vir [M ] 1.4 x 10 15 1.3 x 10 15 1.1 x 10 

At [(xlO 9 ) years] b 8 2 1 



a The centers of the clusters are taken at the electron density 
peaks. 

b Time interval since the epoch when the mass of the cluster 
had half of its final virial mass. 
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by expectations of the primordial or galactic-wind induced initial seed fields, respectively. 
Then, the evolutionary history of the magnetic field was followed under the ideal magneto- 
hydrodynamic approximation. Several interesting results were revealed in their paper. First, 
the initial field configurations are not important for field configurations in the final stage of 
cluster evolution. Secondly, the final field strength is amplified to reach 1/iG by gravita- 
tional collapse and shear flows which induce Kelvin-Helmholtz instabilities. Although the 
simulation ignored effects associated with individual galaxies in the cluster on the evolution 
of the magnetic field, it is worth mentioning that the simulation results can explain some 
observational implications of the rotation measure (also see Dolag et al. 2001a,b). 

In the following, we consider three clusters for demonstrating the magnetic field recon- 
struction, and hereafter we refer to these clusters as model A, B and C. The redshifts z, virial 
radii r vir , virial masses M vir and the interval time since the formation of the clusters, At, are 
given in Table 1. Here, we simply define the formation time as the epoch when the cluster 
had half of its final virial mass. Fig. 1 shows the average electron density (n e ) along the line 
of sight s, (n e ) = 1/(2L) J^ L dsn e , where 2L (= 4 Mpc) denotes the size of the simulation 
volume. The electron density distribution of model A looks spherically symmetric. Since 
the main amplification mechanism in the simulation is gravitational collapse, the magnetic 
field strength should be roughly in proportion to the electron density. Thus, we expect 
that the distribution of magnetic field strength can be also roughly considered as spherically 
symmetric. For this cluster model, therefore, the spherical symmetry approximation, which 
is a basic assumption for the following reconstructions, would be appropriate for both the 
magnetic field strength and the electron density. Model B has two density peaks in the 
cluster inner region, and the ratio of these peaks is ~ 1 . Hence, as a result of our analysis 
which takes the electron density peak as the cluster center, the reconstructed quantities may 
deviate from the spherically averaged true values around the second peak. Model C has a 
subcluster in the outskirts region, and the spherical symmetry assumption may not be ap- 
propriate for this cluster around the subcluster. Thus, the B and C clusters are considered 
to demonstrate that our method can work well for clusters with substructures. Note that the 
redshift of C cluster is higher than those of the other two clusters, which is the one reason 
for the elongated structure. The simulation sequence indeed shows that the subcluster is 
merged with the main cluster component at the lower redshift. It should be again stressed 
that we will take the cluster centers as the electron density peaks in the following analysis. 

The ratio (3 between thermal and magnetic pressure in the simulated clusters is specifi- 
cally addressed in Dolag et al. (2001a), where it is shown in Fig. 1 that the (3 parameter is 
almost constant within the virial radius. 
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3. Rotation measure and primary polarization 

3.1. Small scale limit approximation of the Stokes parameters 

Since the Faraday rotation effect due to the magnetized clusters on the CMB polar- 
ization is important only on the small angular scales, we can safely employ the small angle 
approximation (Zaldarriaga & Seljak 1998; Takada et al. 2001). In this limit, the Q and U 
fields of the Stokes parameters can be expressed using the two-dimensional Fourier transfor- 
mation in terms of the electric (E) and magnetic (B) modes (Kosowsky 1996; Zaldarriaga, 
& Seljak 1997;Kamionkowski, Kosowsky, & Stebbins 1997 a,b; Hu, & White 1997 a,b) as 

^[£ z cos20 z -5 z sin20 z ] e d 6 , 

j^- 2 [E l sm2<j )l + B l cos2<j )l \e d - e , (1) 

where I and 0; are defined as I = Z(cos0 z , sin0 z ) in the fixed (x,y)-orthogonal coordinates. 
Here, E} and B} are the Fourier coefficients for the primary E- and -B-modes, respectively. 
We use the CMBFAST code (Seljak & Zaldarriaga 1996) to generate E^ and ignore the B- 
mode for simplicity, since the vector- and tensor-mode primordial fluctuations that induce 
the S-mode have not so far been detected. Thus, under the assumption of Gaussian random 
fields, we can easily make a primary polarization map (e.g., see a method of making the 
primary CMB temperature map in Takada & Futamase 2001 in more detail). 



3.2. Faraday rotation 

If a linearly polarized monochromatic radiation of frequency v is passing through a hot 
plasma in the presence of magnetic fields along the propagation direction s, the polarization 
vector will be rotated by the angle (the rotation measure), Arm (see e.g., Rybicki & Lightman 
1979); 

ArM = ^c^S^^'^ 

-^^^^(^V^^)^)' (2) 

where Ao is the observed wavelength, and e and m e denote the electron charge and mass, 
respectively. This Faraday rotation effect should rotate the primary CMB polarization vec- 
tor as a secondary effect. From eq. (2), one can see that the rotation measure arises from 
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the combined contributions of the electron density n e and the line-of-sight component of the 
magnetic fields; B s = B ■ s. To break this degeneracy, we first consider a way to recon- 
struct the electron density fields in the next section. It is worth noting that the wavelength 
dependence could be used to separate the Faraday rotation effect from the primary CMB 
polarization and other secondary effects. Fig. 2 shows an example of simulated CMB po- 
larization maps with and without the secondary effect due to the Faraday rotation effect, 
where we have considered cluster A model. 
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Fig. 2. — A example of the simulated map of the primary CMB polarization (white lines) 
and the map (black lines) including the Faraday rotation effect for cluster A model. The 



map is overlapped with the average density fields along the line of sight, (n e (cm 
1/(2L) J^ L dsn e as in Fig. 1. The axes x and y show the angular scale in degrees. 



)) 



4. Reconstruction of n e 

To break the degeneracy between B s and n e in the rotation measure (2), we first consider 
a method for reconstructing electron density fields from possible observed maps of the A-ray 
surface brightness and SZ "flux" , which can be produced directly from the simulation data of 
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the clusters. We expect that typical clusters can be reasonably approximated as spherically 
symmetric bodies. Then, we can readily write down the X-ray surface brightness, Sx, and 
the Compton y parameter for the SZ effect, Sy, at a given frequency band: 



where s = \/r 2 — r' 2 , r' and r denote the projected separation and spatial radius from the 
cluster center, respectively, and T e is the electron temperature. Although we will use the 
physical scale r' throughout, Sx and Sy can be also given by the angular separation 9 from 
the center via the relation r' = cIaO, where cIa is the angular diameter distance to the cluster. 
Although A x and A Y can be computed once the observational frequency bands are specified, 
their explicit values and expressions are not relevant for the following analysis. Here, a(T e ) is 
the X-ray emission coefficient. We focus on the thermal bremsstrahlung only for simplicity 
and set a oc \fT~ e , although in reality one should take into account line emissions as well. 

Under the spherical approximation, we can employ the Abel's integral to reconstruct 3D 
electron density fields following the method developed by Yoshikawa & Suto (1999). Let us 
briefly summarize the concept. As in eq. (3) and (4), observable quantities f(r') for clusters 
are often written as integrals of the corresponding three-dimensional quantities g along the 
line of sight as 



If the quantity g depends on the radius r only (i.e. g — g(r)), we can use the Abel's integral 
to find 



Thus, we can use this transformation formula to deproject the quantities in the clusters, 
assuming spherical symmetry. Applying the transformation above to Sx(r') and Sy(r') 
given by eq. (3) and (4) yields 




(3) 



(4) 




(5) 




(6) 





and 
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From the above two equations, the electron density can be reconstructed as 



(l/(A x n) JdS x (r')r'/Vr' 2 - r 2 ) 



2/3 



n e (r) 



(l/(A Y n) JdSY(r J )r'/Vr' 2 - r 2 ) 



1/3- 



(9) 



In the following, we demonstrate the performance of the above reconstruction technique. 
Since the simulated clusters are not spherically symmetric, we cannot directly apply eq. (9) 
to the simulation data. Therefore we performed the following procedure along the lines of the 
Abel's integral concept. We first computed the circularly averaged values of the simulated 
Sx(r') and Sy(r') in the annulus of a given projected radius r'. We then employed eq. 
(9) to the averaged S x and S Y , giving an estimation of the spherically averaged electron 
density profile. The left panel of Fig. 3 shows the averaged values of Sx(r') and Sy(r') as a 
function of the projected radius r' calculated from the simulated maps of cluster model A, 
where the error denotes the 1 a dispersion from the average value in each annulus (see also 
Yoshikawa & Suto 1999). Note that we do not attempt to include observational errors. The 
right panel of Fig. 3 plots the reconstructed electron density profile with la error bars for 
cluster A. For the error assignment, we have constructed 200 realizations of the S x (r') and 
SV(r') assuming Gaussian distributions of the errors at each bin, and then we calculated 
the average and la dispersion from the reconstructed n e in the realizations. It is apparent 
that the reconstructed density profile well matches the true spherically averaged profile of 
n e directly computed from the simulation data. This is partly because cluster A model has 
no distinct substructure and the spherical approximation can be valid (see Fig. 1). 



Although one may imagine that the Abel's integral can be similarly applied to the 
reconstruction of the magnetic fields, it is not so easy as the reconstruction of the electron 
density fields, since the magnetic field is a vector quantity and various observations have 
implied the existence of reversal scale of the intracluster magnetic fields (e.g., Carilli & 
Taylor 2002). The magnetic fields can be treated like scalar quantities on scales below the 
coherence length where the directions of the magnetic fields are the same. We therefore 
expect that the Abel's integral method can be effectively applied in order to estimate a 
typical field strength for the coherence scale. For this reason, we pay attention to developing 
a method for reconstructing the magnetic field strength as well as a way to estimate the 
coherence length from the observable quantity, i.e., the rotation measure. It should be again 
stressed that the magnetic field strength and the coherence length can be reconstructed as 
a function of radius from the cluster center, which benefits from the fact that the CMB 
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Fig. 3. — In the left panel, the upper line shows the normalized SZ ^-parameter fields, Sy, 
and the lower line denotes the normalized X-ray surface brightness, Sx, for cluster model 
A. The plotted values are the average values of those simulated quantities in the annulus of 
a given projected radius r', while the error bars show the la deviation in the same annulus. 
The horizontal axis denotes the radius r normalized by the virial radius r v ; r . The right panel 
shows the reconstructed electron density n e with ±1<t errors (see text in more detail) and 
the true value of spherically averaged density profile directly obtained from the simulation. 

polarization field is a continuously varying field on the sky, if the secondary effect can be 
measured. 



5.1. The dispersion of the rotation measure 

To apply the Abel's integral to the observed rotation measure, we consider a simplest 
model of tangled magnetic fields; we assume that the magnetic field strength and the coher- 
ence length depend only on the radius from the cluster center. This is an assumption of our 
model, but our results show that it works quite well, even despite the complicated structure 
of the coherence length within the simulated clusters. Even for this case, we should gener- 
ally consider the random orientations of the magnetic fields for each cell with the coherence 
length, whereas the rotation measure along any line of sight can be regarded as a random 
walk process. Hence, the rotation measure in the direction of the projected radius r' can be 
expressed as a sum of the contributions from N cells along the line of sight: 

N 

A RM (r') = y^n e ,Jnl-B w l COSXn- (10) 

n 

where x represents the angle between the magnetic field direction and the line of sight. 
Following the method developed by Lawler & Dennison (1982), the expectation value of the 



rotation measure in the annulus is zero and the 1 a dispersion, a rm {?')■, can be estimated as 



N M 



a RM ( r ') = ^2^2( n e 

\B n \ cosxn x n e \B m \ cos 

Xm/i 



n 



rn 



^ ^ i^n\B n \n etn ) , 



n 



2 f°° rdr 



l(r)n 2 e (r)B 2 (r), 



(11) 



3j r > Vr 2 - r' 2 



where in the second equation we have used (cos 2 Xn) = 1/3, assuming that the magnetic 
field orientations in neighboring spherical shells are uncorrelated. Unordered gas motions 
imprinted during gas formation stir the magnetic fields in the simulated clusters well enough 
for its auto-correlation length to be relatively small, of order 50 — lOOkpc/h. This is in good 
agreement with observations of the Faraday rotation patterns (Kim et al. 1990, Feretti et 
al. 1995). Since we can have a sufficient number of cells at a large separation r', eq. (11) is 
likely to give a good approximation because of the central limit theorem. Interestingly, eq. 
(11) allows us to use the Abel's integral (see eq. (6)) in order to obtain the field strength: 



This equation means that the magnetic field strength can be estimated from the rotation 
measure once the coherence length and the electron density profile are given. Therefore, the 
next thing we should do is to consider a method for estimating the coherence length from 
the rotation measure fields. 



The two-point correlation function will be a useful and simplest quantity to estimate 
the coherence length of the magnetic fields. Let us first consider the two-point correlation 
function of the magnetic fields, which can be directly computed using the simulation data. 
From the meaning of the rotation measure, we here concentrate on the line-of-sight compo- 
nent of the magnetic fields, B s . The the two-point correlation function of B s , £ Bs , can be 
properly defined as 



where x and x + y represent the three-dimensional position vectors from the cluster center 
and (...) denotes the average among all possible pairs within the considered cluster subjected 




(12) 



5.2. Coherence length of magnetic fields 




(13) 



to the conditions of \x\ = r and \y\ = r' . Note that, from our definition of £, we have 
£_B s (r' = 0|r) = 1. Around a radius r from the cluster center, the two points, which are 
separated by more than the coherence length /, are considered to be uncorrelated because of 
the random field orientations, leading to C,b s ~ for r' £ I. Then, the coherence length at a 
given radius r will be roughly estimated from r' that satisfies the condition A^B 3 (r|r') = — 1. 
We can thus define the coherence length of magnetic fields as a form of the differentiation of 
the correlation function: 



Likewise, we can define the two-point correlation function of the rotation measure fields 

as 



where x represents the projection (2D) position vector in a circle of the radius r from the 
cluster center, and y represents a point in a circle of the radius r' from the vector x. Although 
the correlation function of the rotation measures contains some information on the magnetic 
field coherence length, we should bear in mind that the magnetic field in the rotation measure 
is included not as a single field at a fixed point but as a summation of the fields weighted 
by the electron density along the line of sight. Thus, the two correlation functions of the 
magnetic fields and the rotation measures, £b 3 and £_rm", do not necessarily match each other. 

Fig. 4 plots £ Bs and £ RM as a function of r' for cluster A model. We here show the 
results for r/r vir = 0.005, 0.20, 0.30 and 0.55. Although it is not theoretically apparent 
whether these correlation functions match each other, these correlation functions seem to 
share the same properties. In fact, as shown in Fig. 5, the two coherence lengths that are 
estimated using eq. (14) from and £rm have almost the same features. This success can 
be explained as follows. The rotation measure arises mainly from the region near the cluster 
center because both the electron density and the magnetic field strength are relatively large 
there in the simulation data. Hence, the rotation measure is likely to reflect the property of 
the electron density fields and the magnetic fields on the plane which contains the cluster 
center. The magnetic fields determine the sign of the rotation measure because of the 
orientation dependence, and the electron density field is relatively a more smoothly varying 
function with respect to a radius than the magnetic fields. Thus, we could neglect the 
contribution of the density fields to the properties of the two-point correlation function of 




(14) 



Here, to determine l(r) from the simulated clusters, linear fitting is applied where \ < £ < 1 
in practice. 




(15) 
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the rotation measure, and expect that the following relation roughly holds: 



(16) 



In fact, we will see in §5.3 that this relation is also valid for B and C models of clusters. 
Then, without any additional observational quantities, we can reconstruct the magnetic field 
strength using the coherence length derived from the rotation measure, combined with the 
reconstructed electron density fields (see eq. (12)). 



We have presented a method for reconstructing the magnetic fields using the Abel inte- 
gral (12) under a restricted spherically symmetry assumption. In what follows, we demon- 
strate the performance of our method by comparing the reconstructed results with the true 
properties of the magnetic fields directly calculated from the simulation data. From eq. (12), 
we expect that our method gives the spherically averaged profile of |-B(r)| for a given radius, 
which is shown in the upper panel of Fig. 6 for cluster A model as in the left panel of Fig. 
3. One can see that the averaged profile of -B(|r|) is relatively noisy. This implies that the 
spherically symmetric assumption does not accurately hold for the simulated cluster. For 
this reason, we will also consider a following cumulative quantity as a radial profile: 



Note that, to obtain the reconstructed profile of is(r), we first reconstruct B(r) using eq. 
(12) and then compute /b(?") using the equation above. The profile of is(r) for cluster A 
model can be shown in the lower panel of Fig. 6, which implies that Ib is not so different 
from B(r) because of its outer shell weighted form. The quantity Ib{ t ) clarifies the tendency 
of the increase or decrease of the magnetic field strength due to its cumulative form. 

Fig. 7 shows the reconstructed electron density fields, coherence lengths, magnetic field 
strengths and their cumulative quantities (Ib) for cluster models A, B and C. 

Cluster A model is sufficiently virialized and we can consider spherical approximation 
as a reasonable assumption. As expected, all quantities of this cluster are well reconstructed. 
One can also see that the coherence length increases with radius r and the coherence length 
at the cluster center is small due to the large random motion of the fluid. The coherence 
length turns out to be an important property of the magnetic fields for reconstructing the 
magnetic field strength precisely from the rotation measure. 



5.3. 



Results 




f<m drr 2 \B(r)\ 
J Q r drinr 2 



(17) 
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Cluster B model has two density peaks, and thus should be less virialized than the model 
A. Around the second peak (r/r vir ~ 0.15) of this cluster, the electron density is reconstructed 
larger than the true value directly calculated from the simulation data, because the circularly 
averaged X-ray and SZ fields in the annulus including the second density peak tend to give 
a reconstructed electron density larger than the spherically averaged density fields around 
the radius of the second peak (see eq. (9)). The coherence length of the rotation measure 
relatively matches the value of the magnetic field computed from the simulation data. Only 
around the second density peak where the random motion is large, however, the coherence 
length is reconstructed small relative to the true value. This may be because the random 
motion effect in the subcluster decreases the coherence length of the circularly averaged 
rotation measures more effectively than that of the spherically averaged magnetic fields. 
It should be noticed that the increase of the coherence length with radius is weaker than 
in cluster A. This may be due to the fact that in cluster B, which is less virialized than 
cluster A, the magnetic fields are not so efficiently tangled by the shear flow as in cluster A. 
Although the magnetic fields are reconstructed with the density fields which deviates from 
the true value, the magnetic field strength seems to be well reconstructed. This is because 
the rotation measures which include the projected density fields show the same increase as 
the reconstructed density fields around the density peak, and the magnetic field strength is 
constructed from the ratio of these two values, which cancels out the increase of the density 
fields due to the second peak. 

Cluster C model has a subcluster in the outskirts, since the cluster is at the highest 
redshift and dynamically youngest system than the other two cluster models (see Fig. 1 
and Table 1). In this cluster, the subcluster raises the reconstructed electron density near 
f/fvh ~ 0.7 as the result for cluster B around the second peak. The coherence length is 
reconstructed smaller than the true value around the subcluster, because the assumption of 
spherical symmetry is not appropriate as the result for cluster B model. We again stress that 
the relatively small slope of the coherence length with radius may be due to the insufficient 
virialization. In this case, the substructure in the outer region leads to the result that the 
magnetic field strength is reconstructed larger than the true value. This overestimation of 
the magnetic field strength is also caused by the deviation from the spherical symmetry due 
to the concentration of the magnetic field in the substructure. Nevertheless, it is interesting 
that the magnetic field strength is relatively well reconstructed at the inner range of a radius. 

It is worth clarifying why the magnetic field strength in cluster C is relatively stronger 
in the outer region than the other models. The reason was discussed in Dolag et al. (2002): 
The amplification of the relatively weaker magnetic fields in cluster outskirts by adiabatic 
compression and shear flows is more efficient than in cluster cores; firstly because the gas 
flow is more radially ordered outside the core, thus amplification by shear flows is more 
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efficient; secondly because the Alfven speed is lower in weaker fields and thus closer to the 
velocity differences across flows, which also makes shear flow amplification more efficient; 
and thirdly because the stronger fields in cluster cores become dynamically important and 
thus counter-acts further amplification. This explains why magnetic field amplification can 
well be more efficient in dynamically more active, younger systems like cluster C. 

In total, the reconstructed profiles of the strength and coherence length of the magnetic 
fields agree with the true values for the three cluster models, in which the two clusters have 
the substructures in the inner and outer regions, respectively. 

6. Conclusions 

In this paper we have constructed a new method for reconstructing magnetic fields 
in galaxy clusters from the Faraday rotation effect on the CMB polarization, combined 
with possible observed maps of the X-ray emission and SZ effect on the CMB. Our results 
imply that the CMB polarization can be potentially used to reconstruct detailed radial 
profiles of the coherence length and strength of the magnetic fields. It was shown that the 
coherence length estimated from the rotation measure matches that of the magnetic fields. 
Therefore, we do not need any other information than the rotation measure for estimating 
the magnetic field coherence length. This coherence length is not only an important quantity 
for determining the magnetic field strength but also could reveal the origin of the initial seed 
fields. 

To reconstruct the field strength, we consider the dispersion of the rotation measure 
fields in the annulus of a given projected radius, motivated by the random walk process 
caused by random orientations of the magnetic fields for a cell with the coherence length. 
However, the deviation of the field strength from the spherically averaged value may also 
increase the dispersion of the rotation measure, which we do not consider in this paper. This 
increase of the dispersion should lead to the increase of the reconstructed magnetic field 
strength. Nevertheless, the field strengths are well reconstructed in the three cluster models, 
and are reconstructed within a factor of a few even around the area where the dispersion of 
the field strength is large. Anyway, we expect that the method constructed in this paper 
will be a powerful tool for probing the intracluster magnetic fields. 

Finally, we comment on the feasibility of this method. It is a great challenge for current 
technology to detect the secondary effect of the intracluster magnetic fields on the CMB 
polarization. The rotation angle becomes about 1 — 10° at a frequency of 10GHz (e.g. 
Arm ~ 1 — 10°(10GHz/z/) 2 ). The sensitivity of the detector, which is needed to detect 
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the Faraday rotation, is of order 1 /j,K. The angular resolution needed to reconstruct the 
structure of the magnetic fields is estimated from the minimum coherence length of the 
magnetic fields in the simulation cluster as ~ 50kpc ~ 20" (at z ~ 0.1). The frequency 
dependence of the Faraday rotation can be also used to discriminate the effect from other 
secondary signals. Many observations are ongoing and planned for measuring the CMB 
polarization. We expect that future extensive observations of the CMB polarizations will 
allow reconstructions of intracluster magnetic fields with sufficient accuracy, which should 
give a crucial key to understanding the origin of intracluster magnetic fields. 
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Fig. 4. — The upper panel shows ^B s (r'|r) for cluster A model, and the lower panel shows 
^rm(^'I^) for the same cluster. The correlation functions are plotted for each radius r which 
is the distance from the cluster center. The horizontal axis denotes the radius r' normalized 
by the virial radius r vir . Note that the radius parameters r' for £ Bs and £rm are distances 
from a fixed point located in a sphere or circle of the radius r from the center of the cluster, 
respectively. 
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Fig. 5. — The coherence lengths of the magnetic fields in cluster A. The horizontal axis 
denotes the radius r normalized by the virial radius r vir . The dotted and solid curves show 
the results estimated from £g s and using eq. (14), respectively. Those coherence lengths 
well match each other within the considered range of radius r. 
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Fig. 6. — The upper panel shows the spherically averaged profile of \B(r) | for cluster A model, 
where the error bars denote la deviations from the average value for each bin of radius r 
as shown in Fig. 3. Similarly, the lower panel shows the cumulative field strength I B (r)[G], 
which is computed using eq. (17). The horizontal axes denote the radius r normalized by 



the virial radius r vir . 
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Fig. 7. — These panels show the comparisons of the reconstructed profiles of the electron den- 
sities, the coherence lengths, the magnetic field strengths and the cumulative field strengths 
with the true profiles computed from the simulation data. The left, middle and right panels 
are the results for clusters A, B, and C, respectively. Those quantities are reconstructed 
using equations (9), (14), (12) and (17), respectively. 



